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ABSTRACT. 

We calculate structure functions at small x both under the assumption of a 
hard singularity (essentially, a power behaviour x~ x , X positive, for x — > 0) or that of 
a soft-Pomeron dominated behaviour, also called double scaling limit, for the singlet 
component. A full next to leading order (NLO) analysis is carried for the functions 
F2, fbiue and the longitudinal one Fl in ep scattering, and for £^3 in neutrino scatter- 
ing. The results of the calculations are compared with experimental data, particularly 
the recent ones from HERA, in the range x < 0.032, 10 GeV 2 < Q 2 < 1 500 GeV 2 . 
We get reasonable fits, with a chi-squared/d.o.f. around two unities, with only three- 
four parameters for both assumptions, but none of the assumptions is by itself able to 
give a fully satisfactory description of the data. The results improve substantially if 
combining a soft and a hard component; in this case it is even possible to extend the 
analysis, phenomenologically, to small values of Q 2 , 0.31 GeV 2 < Q 2 < 8.5 GeV 2 , and 
in the x range 6 x 10 -6 Si x 5^ 0.04, with the same hard plus soft Pomeron hypothesis 
by assuming a saturating expression for the strong coupling, 

& 3 (Q 2 ) = An/p log[(Q 2 + A 2 B )/A 2 S ]. 

The description for low Q 2 implies self-consistent values for the parameters in the 
exponents of x both for singlet and nonsinglet components. One has to have, for the 
Regge intercepts, a p (0) = 0.48 and ctp(0) = 1.470 [A = 0.470], in uncanny agreement 
with other determinations of these parameters, and in particular the results of the 
large Q 2 fits. The fit to data is so good that we may look (at large Q 2 ) for signals of 
a "triple Pomeron" vertex, for which some evidence is found. 

The quality of the calculations of F2 , and of the predictions for Fbiuc , Fl is 
only marred by the very large size of the NLO corrections for the singlet part of F2. 
This, in particular, forbids a truly reliable determination of the QCD parameter, A. 



* This paper includes the results from FTUAM 96-39 [hep-ph 96 10380] and FTUAM 96-44 [hep-ph 96 12469] 
'* e-mail: fjy@delta.ft.uam.es 



-Adel, Barreiro and Yndurain- 



1. INTRODUCTION 



In two recent papers' 1 ' 2 ! it was shown how the recent HERA data' 3 ' 4 ! on electroproduction at small x could 
be very well fitted by the formulas, proposed long ago by C. Lopez and one of us' 5 ' 6 ' which to leading order (LO) read, 

F 2 (x,Q 2 ) = (e 2 q )[F s (x,Q 2 )+F NS (x,Q 2 ) 
F S {x,Q 2 ) 



x^O 



B s a- d +x- x , 



F NS (x,Q 2 ) ~B NS a- d » s xe,p~0.5. 

x^O 



(1.1a) 
(1.1b) 
(1.2) 



Here d + and d^s are known quantities related to the singlet anomalous dimension matrix D and to the nonsinglet 
anomalous dimension; Bs, -Bjvs, A are free parameters. In ref. 2, the analysis was extended to the gluon structure 
function, while the longitudinal one (to LO) had been considered in ref. 8. One has, 

U( i .( 

x^O 

and 

a s G'fo: s f 1 t ATpUf Bq 



F G (x,Q 2 ) ~B G a- d +x 



(1.3a) 



R = Fl/F 2 ~ r - 

x^O 7T 



1 + 



(1.3b) 



tt(2 + X){ (3 + X)C F B s 

with Bq/ Bs = (d + — Dn)/D\2- These formulas follow at leading order from perturbative QCD, plus the assumption 
that the leading singularity in n of the matrix elements (for e.g. the singlet case) 



(p\a n \p): a„ = 



q f jd...dq f 
G8...8G 



V 



ra-2 



are located to the right of those of the corresponding anomalous dimensions, D, dws (For more details, see refs. 1, 5, 
6,7). 

An alternate possibility, that may be called "soft-Pomeron" dominated, occurs when the singularity of D is the 
leading one; it was proposed first by De Rujula et a I J 9 ' The ensuing LO behaviour for the structure functions Fs, Fq 
was evaluated in detail by F. Martin' 10 ! an d f or R j n re f 8. The next to leading order (NLO) corrections to Fs were 
given in ref. 11. To LO one has, 



F s ~ 



CO 



log a; I 



9|loga:|logK(g§)/ as (Q 2 



47r 2 (33-2n/) 



xexp< Jd \\ogx\ 



log 



a s (Q 2 ) 



- di log 



a s (Q 2 ) 



(1.4a) 



F G ~ 



9co 



33- 2nt 



5767r 2 |logx|log[ as (g 2 )/a s (Q 2 )] 



x exp < J d 1 logx| 



log 



d a = 144/(33 -2n f ), d x 



a s {Q 2 ) 

(33 + 2n//9)/(33- 2n } ), and 

(33 - 2n f )\ \ogx\ 



- dx log 



a s (Q 2 ) 



(1.4b) 



a s (Q 2 ) 
2ir ' 



(1.4c) 



_log[ as (Q 2 )/a s (g 2 )]_ 

to be compared with Eqs. (1.1), (1.3). We will discuss this possibility in Sect. 4. 

Returning to our first case, the analysis was extended in ref. 2 to practically the whole range of HERA data at 
the cost of introducing phenomenological correction terms for "large" values of x, x > 0.01, and small Q 2 < 12 GeV 2 . 
Moreover, only the LO prediction for R was evaluated. Here we go a few steps forward, in the following directions. 
First, we perform a full NLO analysis, including that of the longitudinal structure function. Second, we extend the 
analysis to incorporate the function xF?, for neutrino scattering which provides pure nonsinglet function, hence a 
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combination independent from that in (1.1a): this helps stabilize the results. Then, for the hard Pomeron case, we 
include theoretically justified corrections, and resummations, which enable us to extend the range of validity of the 
formulas (as determined e.g. in ref. 1) and to paliate somewhat the effects of the hughe size of the NLO corrections 
to the singlet component of F 2 , which are of some ~ 6a s . 

For the soft Pomeron dominance hypothesis we again perform a full NLO calculation of F$, Fq and Fl. NLO 
corrections are, also in this case, very large; their size is indeed the only serious drawback for our results, otherwise 
able to describe reasonably well the HERA data in a wide range, and providing believable predictions for the gluon and 
longitudinal structure functions. This agreement of theoretical predictions and data is essentially true both for the soft 
and hard Pomeron hypotheses, an apparently surprising fact that is discussed in Sect. 5. Here we argue that a possible 
reason is that Fs contains both a hard and a soft piece. In fact, we are able to give excellent fits to all HERA data 
(x < 0.032) by using a formula sum of (1.1) and (1.4). This provides us with the best fits to the data, the parameters 
of which are reported in Table X. 

In Sect. 6 we show that the fit with a soft plus a hard Pomeron may be phenomenologically extended to low Q 2 , 
down to 0.31 GeV 2 , provided we make a saturation assumption for the strong coupling a s and satisfy a self-consistency 
condition for the parameters A, p. This fixes these parameters to the values A = 0.470, p = 0.522, in uncanny 
agreement with other (in particular high Q 2 , Sect. 5) determinations. 

2. THEORETICAL EVALUATIONS (HARD SINGULARITY) 

We will here briefly rederive the extension to NLO of the equations governing the behaviour of structure 
functions as x — > 0; not only for ease of reference, but because the large size of the singlet NLO corrections makes 
it convenient to use formulas more precise than those employed in refs. 2, 6. We will also extend the analysis to the 
longitudinal structure function at NLO. 

2.1. Nonsinglet 
Defining the moments 

PNs(n,Q 2 )= [ dxx n - 2 F NS (x,Q 2 ), (2.1) 
Jo 

they satisfy the QCD evolution equations 

PNs(n,Q ) = e Jt « C N s{n,a s {Q 1 )), (2.2) 

where t = log Q 2 , g is the coupling constant, 7_/vs the nonsinglet (NS) anomalous dimension, and C/vs the NS Wilson 
coefficient 1 . The first singularity of 7^5 to LO and NLO lies at n = 0. If we assume this to occur to all orders, and 
that Fns behaves like a power x p as x goes to zero, we get the following behaviour from (2.2): 

F NS (x,Q 2 ) ~ (Const.)^e"^o dt ' 7Ns(1 " P ' s(t ' )) C JVS (l-p,a s (Q 2 )). (2.2) 

x^0 

Expanding ~fNS, 9, G'ns to second order and integrating we get, 

F NS (x,Q 2 ) ~ B NS xP 

x^0 

x exp [d NS (l - p) WW 2 ) + gA " (1 ^ WQ2) +-..}, (2-4) 

Bns a constant and 

d NS (n) = -7#k»)/ 2 #) 

Pid NS (n) 7^(n) 
qNs(n) ^^o~ + ^o- 



1 The anomalous dimensions and coefficients are collected in the Appendix for ease of reference. 



1 | C^ s (l-p)a s (Q 2 ) | 

47T 



- 2 - 



-Adel, Barreiro and Yndurain- 



The values of the quantities /3j, Cjvs, may be found in refs. 6, 7. Eq. (2.4) may be conveniently rewritten 
(suppressing the dots) as 



x^a I 47r j 



(2.5) 



a gws(l-p)a s /47r 



)2\] - d Ns(l-p) p 



xe"" 1 ' r [a s {Q 2 )\ '• " x r 

an expression 2 in which the modifications of (1.2) due to the NLO corrections is apparent 
One can also expand the exponent in (2.5) and get 

«ats(1 - p)a s {Q 2 



F NS (x,Q 2 ) ~ B NS {1 + 

x^O 



47T 



a; dNS x p , 



vns = C ( xl + q N s- 

For p = 0.5, a value that follows from a Regge analysis and that we will adopt here, one finds 

VNs\p=0.5.n f =4 = 3.42, 

so the NLO correction is small and we may use (2.5) or (2.6) indifferently. 
2.2. Singlet 

Eq. (2.2) is now replaced by the coupled equations 

fi(n,Q 2 ) = (Const.)T e" ^ " l{n ' 9[t ' )) C{n,a s {Q 2 )) 



(2.6) 



(2.7) 



(2.8) 



fi(n,Q 2 ) 



Vs{n,Q 2 ) = $l&xx n - 2 F s {x,Q 2 ) 
^ G {n,Q 2 ) = fidxx n - 2 F G (x,Q 2 ). 
Here 7, C are square matrices; the operation T in (2.8) is like the familiar time ordering operator but it now orders 
in t = logQ 2 . This ordering, and the matrix character of the equations complicates the singlet analysis, the details of 
which may be found in refs. 6, 7. To next to leading order we may easily write the analogue of (2.4), (2.5) as 



F s (x,Q 2 ) ~ B s x- 

x—>0 



1 | c s (l + X)a s (Q 2 ) | 

47T 



x cxp [d + (l + A) W(Q 2 ) + gg(1+ ^ (Q2) + 



} 



or, suppressing the dots and in a form easier to compare with the LO expression (1.1b), 

x^O 47T J 



xe 



9s (l+A)a s /47r 



[a s (Q 2 )]- a+{1+X) x~\ 



(2.9) 



(2.10a) 



A corresponding equation for the gluon component we will consider later. In above equations we have 



D = - 



r (o) 
2/V 



_ r (i) d+ - Dn (1) 



0id+ , 711 
qs = — — h 



721 



D 



12 



0o 2/3 2p (d- - d+ + 1) d_ - d+ ' 
7 = S- 1 7 (1) S, 



and S is the matrix that diagonalizes D: 



S _1 DS 



d+ 
d- 



Eq. (2.5) corrects the sign misprints in Eq. (2.18) of ref. 6. 
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One can also expand the exponential in (2.10a) and get 

Fs(x, Q 2 ) ^ B s [l + Ws(1+ 4 A ?r K(Q2) } a-* +X ~\ (2.10'a) 

W S = C S + 9s- The 7, g, c, w are collected in the Appendix. 

Unfortunately, is very large. For A = 0.35, Uf — 4, wg — 77.8; for A = 0.47, W5 = 56.7. Therefore, 
we are faced with the choice of using the exponential form (2.10) or the expanded one (2.10'). The exponential form 
has errors of order a 2 because the noncommutativity of 7^, 7^ makes the T-exponential different from the ordinary 
exponential. If we use the expanded form (2.10') we have other errors (also of order a 2 ) due to the large size of 
the neglected term Q[qg{l + X)ot s /Ait\ 2 . It is unclear a priori which of the two procedures will be more accurate, 
although the abnormally large size of qs(l + A) suggests that the exponentiated form will be more precise; note that 
the perturbative expansion still makes sense, for the exponent, in that for reasonably large Q 2 one has 



47T 



< d + log a s 1 . 



In fact, and as we will see, the exponentiated form produces somewhat more satisfactory results than the expanded one. 
At any rate, we will use both (2.10) and (2.10'): one may take the difference as an indication of the theoretical error 
of our calculation. 

Similar considerations of course apply to the gluon component that we discuss next, although in this case the 
correction is much smaller (~ 15a s /47r) so use of exponentiated or expanded form is essentially equivalent here and 
only consistency with the quark component will make us use one or the other. We then have, 

FcM ^ B A^(i' CG(1+A)a ' ( ° !) 



where now 



x^o D 12 { 4tt J (2.10b) 

xe «,(i+A)a./4* [ as(g2) ] -d+d+X) a ^ 

M 1 ) , /^(l) ( D 22 - Dll . D 2 1 \ 

_ Pid + 7n_ 72i Dl2 d ~ Z Dl1 

q °~ f3 2(3 20o(d--d+ + l) d--d+ d+-D u ' 



In expanded form, 



(2.10'b) 



F G {x,Q) x% B s — |1 + |MQ)] x , 

w G (l + A) ee c G (l + A) + q G (l + A). 

2.3. The longitudinal structure function 

We normalize the longitudinal structure function Fl in such a way that one has 

fS + ^NS — *L 



It is also convenient to define the quantity R' by 



For x — > the contribution of Fjvs is negligible with respect to that of F$ and we will accordingly neglect it; the effect 
of taking it into account, to LO, may be found in ref. 8. The function Fl may be evaluated in terms of F$, F G . One 
has, 

F L (x,Q 2 ) = J\y \c L s ( Vl Q 2 )F s Q 2 ) + C G (y, Q 2 )F G Q 2 ) } , (2.12a) 
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where the kernels C L are, 



C|(s, Q 2 ) = Ck s (x, Q 2 ) + C^ s (x, Q 2 



C^ s {x 1 Q 2 ) = [AC F x] 



u s {Q 2 ) , mi 



47T 



-NS 



(X) 



a s (Q 2 ) 
4ir 



+ . . 



c L PS (x,Q>) = c p s L (x) (sagl 



(2.12b) 



C^(x,Q 2 ) = [16n f T F x(l - x)} 



a s (Q 2 ) , (i)L, , (a s {Q 2 ) 



47T 



47T 



+ ... 



The functions (^nsps G are described in the Appendix 13 



Under our assumptions, Eqs. (2.10), the behaviour of Fl follows immediately; we have, 

F L (x,Q 2 ) ~ F s (x,Q 2 ) [ dxx x {C L s { Xl Q 2 ) + V (Q 2 )C^(x,Q 2 )} , 

Jo 



(2.13) 



and T] is the ratio 



V (Q 2 )=lun FG ^ 



+oF s (x,Q*) 

To LO a simple evaluation gives! 8 ! R(x,Q 2 ) ~ i?(°)(x,Q 2 ) with 

i2<°>(s,Q 2 ) ~ r (l + A)^, 



ro(l + A) 



1 + 



47>n/ d 



+ (l + A)-£>n(l + A) ] 
Di 2 (l + A) J" 



(2.14) 



2 + A t (3 + A)C F 

To NLO the calculation is made numerically. For this, define the integrals (whose values may be found in the Appendix), 

dxx x 4 ] s L (x) = C F (C A - 2C F )I S i(X) + C 2 F I S2 (X) + C F T F n f I ss (X), 

pi 



Jo 



f 

Jo 



J 

Jo 



dxx A c y p'g{x) = C F T F n f I PS {\), 



(2.15) 



dxx 



(x) = n f T F [C F I G1 (X) + C A I G2 {\)], 



and write, in exponentiated form, 
Then, 

with J?(°) as above and 

C F a s (Q 2 ) 



*o 2 + A 7T 



R' = + RfW 



[C A Isx + C F (I S 2 - 2I S i) + n f T F (I S3 + I PS )} g 



+ 



ArifTp 
(3 + \)C F 



1+ (2 + A)(3 + A) (CrJoi+c ^ 
lb 47t 



(2.16a) 
(2.16b) 

(2.16c) 



3. NUMERICAL RESULTS (HARD SINGULARITY ONLY). 
3.1. The function F 2 . 

3 These quantities were first evaluated in refs. 12, 13. Correct values, checked at least in two independent calculations, 
are given in Eqs. (8), (9) of ref. 12 for c^g PS and in Eq. (9) of ref. 13 for Cq\ The value of this quantity given in 
ref. 12 contains an error. The first moments of the c may be found in ref. 14; they are useful, among other things, 
to check the integrals (2.15) here. 
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LO calculations. For ease of comparison between LO and NLO evaluations we repeat here the results of a 
fit to the old (1993) Zeus data! 3 ', as performed in ref. 1. The calculation is carried for 32 points in the range 

x < 1(T 2 , 12 < Q 2 < 90 GeV 2 . 

Because of the size of the experimental errors a LO calculation is sufficient, and the NS contribution may be neglected. 
The QCD parameter A is fixed to 0.2 GeV so that a s (m 2 ) = 0.32. The results are summarized in Table I for 
Hf = 4 flavours. The corresponding values of Be, ro are also given. The agreement of the value of A with the figure 
A = 0.36 ± 0.07 obtained in ref. 3 from data with x > 0.02, Q 2 < 22 GeV 2 is noteworthy. 



Table I. "Old" Zeus data. LO. n f = 4; A(l loop, n s = 4) = 0.200 GeV; ot s {m 2 T ) = 0.32. 



d+ 



(e 2 q )B.s 



Bg/Bs 



ro 



X 2 /d.o.f. 



0.38 ±0.01 2.41 ±0.1 (2.70 ± 0.22) x 10 -3 20.56 ± 0.54 6.24 ± 0.24 ^| 



NLO evaluation. If we only fit the HI points^ with x < 10 2 , Q 2 > 12 GeV 2 using the exponentiated 
formulas, we get a x 2 /d.o.f. is less than one, with parameters reported in Table II. 



Table II. HI plus v data, x < .01, n f = 4; A(n s = 4, 2 loop) fixed to 0.11 GeV. 



(e 2 q )B s {e 2 q )B NS X 2 /d.ol. 



0.3218 1.423 x 10" 



0.390 



However, it is still not possible to give any value for the QCD parameter A. The reason is that the interplay 
between singlet and nonsinglet parts compensates the effect of varying A. For example, a x 2 /d.o.f. less than one is 
attained for 3 MeV < A < 260 MeV. This is why we do not give errors in the parameters in Table II. 



log Q (GeV 2 ) 



7.5 



We may improve the situation as follows. First, 
and as discussed in ref. 1, we can include more points 
limited by a certain Q 2 (x) beyond which corrections to 
the leading behaviour become important. To be precise, 
we choose the HI points with (Fig. 1) 

Q 2 < 150 GeV 2 , for x = 0.013, 
Q 2 < 90 GeV 2 , for x = 0.02, (3.1) 
Q 2 <60GeV 2 , for a; = 0.032, 

with a total of 77 points. Secondly, we incorporate small 
x data (a total of 10 points) from the neutrino struc- 
ture function^ 15 ! xF% which is pure nonsinglet and hence 
provides the independent measurement necessary to dis- 
entangle the singlet and nonsinglet components of F^. 
this, as we will see, gives stability to the results. 

The outcome of the fits is given in Tables III, IV 
with A a free parameter. The x 2 /d.o.f. is reasonable, 
although its increase beyond unity reflects the fact that 
the subleading effects are substantial for the points x = 
0.013 - 0.032. 



Figure 1. Area described by (3.1), bounded 
by a continuous line. Broken line: extended 
region for fit with (3.2). 
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Table III. HI plus v data; x given by (3.1). rif — 4, two loops. 



A 



(e 2 g )Bs 



(e 2 q )B NS xVd-O.f. 



O.O8OIH35 GeV 0.3243 ± 0.0065 1.32llg;^? x 10~ 4 0.254±°;°^ §^ 



Table IV. HI plus v data; x given by (3.1). n/ = 5, two loops. 



A(nf = 4, 2 loops) 



0.1lt° : jg GeV 



(e 2 g )Bs 



(e 2 q )B NS xVd-o.f. 



0.331 ±0.006 1.118lg;SS x 10- 4 0.311±g°? 4 ±22_o 



From these results it is clear that the data do not discriminate between n/ = 4, 5, although the first value is 
slightly favoured. For this reason we will give almost exclusively fits with n/ = 4. 

The values of A we obtain are compatible with standard ones^ 16 !, albeit on the small side. Because the 
parameters are very strongly correlated the errors given are obtained not by varying the parameters independently, but 
by varying only A and treating the other parameters as dependent quantities. It is also important to realize that the 
errors in Tables III, IV and indeed in practically all the evaluations, are purely nominal in the sense that we have not 
taken into account theoretical errors, which are much larger. In fact, the central values for the parameters, especially 
Bs, A, depend very strongly on the theoretical assumptions 4 made; for example, they vary way beyond the nominal 
errors from LO to NLO: compare e.g. Table I with Table II. 

The results reported in Tables III, IV were obtained with the exponentiated formula. If we use the expanded 
one, Eq. (2.10'a) we find the results of Table IV . 



Table IV'. HI plus v data; x given by (3.1). n f = 4; A(n f = 4, 2 loops) = 0.080 GeV. 



(e 2 q )B s (e 2 q )B NS x7d.o.f. 



0.3408 2.54 x 10~ 4 0.268 



99.5 

87-4 



We do not give errors, but the results of the fit for a representative value of A: that for which the fit with 
the exponentiated formula is optimum. This is because there is no optimum reasonable value of A if using the 
nonexponentiated expression; the \ 2 decreases slowly with A down to a few MeV. 

Corrections. Let us now turn to the corrections that will enable us to extend the calculation to all points 
with x < 0.032, Q 2 > 12 GeV 2 (Fig. 1). We take them into account semi-phenomenologically 5 by replacing (2.10a) 
with 

cscts 



JW.,<JVa.{i + 2*} 



(1-x) 



HQ') 



(3.2a) 



and fixing v{Q 2 ) so that, for small Q 2 , we agree with the result of the counting rules for x 
satisfy the momentum sum rule, 



dxF s (x,Q 2 ) 



3n f 



Q 2 ^oo 3ny + 16 



1 and, for large Q 2 , we 



4 This is in fact the reason why we have given a variety of evaluations, and not just the best ones: to get a flavour 
for the systematic theoretical uncertainties. 

5 The similitudes and differences with the more phenomenological procedure of ref. 2 should be apparent. 
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Specifically, we choose 



HQ 2 ) = + 



/ B s [l + c s a s /47r]e^^/ 4 "r(l + A)[16 + 3n f ] \ 



1/(1-A) 



3nja 
vq = 7. 



(3.2b) 



Note that this does not introduce any new parameter. 
For Fns we replace (2.6) by 



'ATS — &NS 



B NS 1 + 



VNSUs 

Air 



}aJ dNS x p (l - x) UNS , 



(3.2c) 



but, because the NS component is only relevant at small values of Q 2 we fix Vns = 3 independent of Q 2 (actually, 
the x 2 varies by less than one unit for < vns < 4). Then, we still write F^ = (e 2 )[Fs + Fns], and, for neutrino 
scattering, 



(3. 2d) 



= 0.00032 



x = 0.0005 



10 12 14 16 





10 15 20 25 rf 10 18 26 34 rf 10 22 34 46 rf 






10 28 46 64 rf ' 10 60 110 160 rf ' 10 90 170 250 ef- 




ts 110 220 330 rf 



x = 0.013 




180 360 540 rf 



x = 0.02 



fl' I 



x = 0.032 



0.5, 



250 500 750 rf 



a 330 660 990 



Figure 2a. Comparison of predictions from Eqs. (3.2a, c), Table V 
(A = 0.14 GeV), with HI ep data' 41 for F 2 . 

The results of the fit are presented in Table V, for the exponentiated expression, Eq. (3.2a). There is 
unfortunately no minimum as a function of A: the \ 2 decreases slowly with A. We thus give results only for two 
representative values of this parameter. The pictorial representation of the fit is given in Fig. 2a (for ep) and Fig. 2b 
for neutrino scattering, both for A = 0.14 GeV. 
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Table V. HI plus v data. n f = 4; x < 0.032, 12 GeV 2 < Q 2 < 1 200 GeV 2 . 



A = 0.10 GeV : 
A = 0.20 GeV : 



0.3183 
0.3286 



1.292 x 10" 
2.257 x 10" 



(e 2 q )B NS xVd-o.f. 

0.328 
0.371 



110-3 
141.3 



= 0.015 



= 0.045 



The x 2 /d.o.f. is slightly larger than one. 
Part of the discrepancy is due to the data, some of 
which is clearly incompatible with the rest. Also, 
one may substantially improve the \ 2 if introduc- 
ing a free parameter in the definition of v{Q 2 ), as 
shown e.g. in ref. 2. However, part of the dis- 
agreement is certainly due to rigidity of the the- 
oretical formulas, and to true deviation from the 
model which occur for "large" values of x. We will 
discuss this further in connection with the analysis 
of the Zeus data, and in Sect. 5. 

Finally, the fact that the x 2 decreases with A past reasonable values is an indication that we are getting here 
an effective value for this parameter, which compensates for the large size of the NLO corrections. 

We may consider fitting with the expanded version of the formula for Fs, i.e., with 



Figure 2b. Comparison of predictions from Eq. 
(3.2c), Table V (A = 0.14 GeV), with neutrino 
xF-i data' 151 . 



Fs(x, Q 2 ) ~ -B s jl + iCS+ 4 f as ) a- d + x~\l xy- 



(Q ) 



and now 



fn.c.(Q 2 ) =V0 + 



B S [1 + (c s + <?sW4tt] r(l + A)[16 + 3n f ] 
v = 7. 



1/(1-A) 



The fit deteriorates clearly; the x 2 is now of some 140, for 110-3 d.o.f. and A = 0.10. This shows that the 
exponentiated version of the formulas is to be preferred, as it probably sums at least part of the large NLO corrections. 
Because of this, we will henceforth use only the exponentiated version of the equations. 

We next consider fits to the more recent Zeus data' 4 !. We will make two choices: first, we fit the neutrino 
data, and all ep points with x < 0.01 using the formula (2.10a). The results are given in Table VI. The chi-squared is 
reasonable, as is the value of A. The values of all parameters are compatible with those found from the fits to the HI 
data. The second possibility is to extend the range to x < 0.025 and use Eq. (3.2), fixing A = 0.135. The results of 
the fit are shown in Fig. 3. We do not show the fit to the neutrino data, which does not differ substantially from that 
of Fig. 2b. The x 2 /d.o.f. is now of 226.1/(120-3). This, as the x 2 /d.o.f. reported for the fit of data with x < 0.01 
in Table VI, are larger than their counterparts for the HI data. A glance to Fig. 3 shows that part of the reason is the 
presence in the Zeus data of fluctuations. These are probably due to systematic errors not taken into account in the 
experimental analysis; they become important for very large Q 2 . Thus, and although the Zeus data appear more precise 
than the HI ones for the lower Q 2 range, 6 the last one are more reliable at large Q 2 . Nevertheless, and as noted in 
the comments to the fit to HI data, it is also clear that the theoretical predictions present systematic deviations from 
experiment, very likely due to the extension of the first beyond their range of validity by use of a semiphenomenological 
expression which is not sufficiently flexible; see Sect. 5 for more discussion. Apart from this, the results are good and 
the parameters of the fits reasonably compatible. The value of A is closer to the accepted one. 



This is probably the reason why the value of A deduced from the HI data is less realistic than that obtained fitting 
the set of Zeus. 
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Figure 3. Comparison of predictions for F2 from Eqs. (3.2a, c), A = 0.135, 
with Zeus ep data 



Table VI. Zeus plus v data. n f = 4; x < 0.01. 



A A (e 2 q )B s (e 2 q )B NS x 2 /d.o.f. 

0.135lg;^GeV 0.301 ± 0.025 1.250^;^ x 10" 4 0.3138 ± 0.007 §§^f 

As a final check on the reliability and consistency of the fits we have fitted Zeus data with x < 0.01, not 
including the neutrino data. We get no definite minimum for A, only constrained by A < 0.2 GeV; but we obtain 
values of the remaining parameters compatible with those obtained including xF^, in particular a very reasonable value 

for B NS : 

(e 2 q )B s = 1.8 x 10- 4 , (e 2 q )B NS = 0.35, A = 0.329. 

In this sense we may say that our analysis is sufficiently precise to predict the NS structure functions from F2 only, and 
this in spite of the relative smallness of fjvs- It is however clear that, as already mentioned several times, systematic 
deviations occur, especially large for x ^ 0.01 (cf. Sect. 5). 

3.2. The gluon structure function. 

We give here the parametrizations to NLO for the gluon structure function that follow from our determination 
of the parameters in the previous subsection, for the full set of points corresponding to a set of parameters intermediate 
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between those given in Table V for F2, 



F G {x,Q 2 )=B G 



-0.25a s 



4?T 



e 16.6a 3 /4. as(Q2) 



2\-3. 182-0. 316 



(1 - x) VG 



vg(Q 2 ) 



MQ 2 ) 



f B G [1 - 0.25a ;i /47r]e 16 - 6a °/ 47r J r(l + 0.316)[16 + 3ny] ) ^ 



:U(i 



V 



16a? 182 



J 



(3.5) 



B G = 25.2 x (e 2 q )B s , (e 2 q )B s - 1.226 x 10" 4 , 
NLO, x < 3.2 x 1CT 2 , 12 GcV 2 < Q 2 < 1 200 GeV 2 ; 

n f =4; 4(2 loop, n/ = 4) = 0.14 GeV. 



The corresponding graphs are shown in Fig. 4, where we give both LO and NLO predictions, the LO calculation with 
values of parameters from Table I. 

There are unfortunately no direct measurements of F G with which to compare our calculations. Indirect 
estimates were made by the HI collaboration, by fitting F2, F G with an exact coupled QCD evolution. The comparison 
with our calculations to LO may be found in ref. 1; the agreement is reasonable, and indeed our estimates are more 
precise than the DGLAP calculation, afflicted by large extrapolation errors. 

More information on F G is obtained from the cross-section jp — > J/ipp, to be discussed in Sect. 5.3. 



10 12 14 16 g2 



10 15 20 25 



10 18 26 34 



10 22 34 46 



x = 0.002 




10 28 46 64 



x = 0.0032 



10 60 110 160 



10 90 170 250 



x = 0.008 



110 220 330 




Figure 4. The gluon structure function Fg to LO (broken line), 
Eq. (3.3), and the optimum NLO one (continuous line), Eq. (3.5). 



3.3. Predictions for the longitudinal structure function. 
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NLO, 0{a s ) corrections to the longitudinal 
structure function are unfortunately very large; not be- 
cause of the direct corrections, but due to corrections 
generated indirectly via the large NLO corrections to 
i<2. Indeed, the value of R' is reduced by more than 
a half from LO to NLO. We give in Fig. 5 a plot of 
LO and NLO calculations. Using Eq. (2.14), and the 
parameters A = 0.20, A = 0.38 (Table la) we get 
the LO result, J?(°); and with Eqs. (2.15), (2.16) and 
the figures A = 0.10, 0.20 and A = 0.324 (cf. Ta- 
ble IVa) we find R' (NLO). Also depicted are a few 
representative data. Note that the dependence of the 
NLO value of R on A is very slight, due to cancella- 
tion of various effects. Thus, the lines corresponding 
to A = 0.10, 0.20 in Fig. 5 fall almost one on top of 
the other. 

To get a further indication on the meaning of the results, we have also calculated R 1 from the effective fit 
at low energy of ref. 1, with A = 0.324, and a s to two loops, but without other NLO corrections, taking n/ = 3 
below Q 2 = 12 GcV 2 , and rif = 4 above. This is the intermediate, dashed curve in Fig. 5. Clearly, one would expect 
that the real R' would somehow interpolate between this, at low momentum, and the full NLO curve, for very large 
Q 2 . The predictions should be checked against experiment when, and if, measurements independent of those of F 2 are 
performed at HERA. We have given the predictions for x = 0; the figures would not change much provided x < 10 -2 . 

4 THE SOFT-POMERON DOMINATED MODEL 

As remarked in the Introduction, the results derived in the previous sections assume that the singlet structure 
functions are dominated, at small x, by the singularities of the matrix elements of the quark and gluon operators. We 
may instead hypothesize that these singularities lie to the left of n = 1, and then the small x behaviour is controlled 
by the singularities of the Wilson coefficients. Specifically, this occurs if one assumes that, for all Q 2 below a certain 
Qq of the order of a typical hadronic scale (say ~ GeV), cross sections behave according to a standard soft-Pomeron 
dominance, 

0tot(Qo, s) ~ Constant, 

for Q 2 < 1 GcV 2 . 

We can then get the structure functions for small x by evolving with QCD the expressions corresponding to 

this,! 9 - 10 - 11 ! 

Fi{x,Ql) ~ a, i = S,G. (4.1) 

We will present a sketchy derivation of the resulting formulas, to NLO. This is of interest because we use the moments 
method, instead of the Altarelli-Parisi equations employed in ref. 11, so we have a nontrivial check of the calculation 
there. 



R(0,Q 2 ) 




80 GeV 2 160 GeV 2 240 GeV 2 



Figure 5. Predictions for R'(0,Q 2 ) to LO 
and NLO. Hatched box: preliminary result 
from Hi' 17 '. Dots: data from ref. 18 (actu- 
ally, at x ~ 0.05). Discontinuous line: inter- 
mediate calculation, see main text. 



4.1. Theoretical calculations: F$ and Fq 

Our starting point is the following relation, proved in ref. 6 to NLO, 

a s (Q 2 ) 6 (l + ar)S _1 (l - flC (1) )/i(n, Q 2 ) = b = independent of Q 2 . 



(4.2) 
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Here a = a s (Q 2 )/4:iT, is the matrix of NLO corrections to the Wilson coefficients 7 and 







a s 
a s 



1 

Wo 



712 (n) 



A 



^n(n) + 2/3ld+ (n) ^ +(n) _ rf _ (n) + 1 
L(n)-t(n) + l %"(») + 



(4.3) 



with S, 7 defined before; explicit expressions for these quantities may be found in the Appendix. Here we only give the 
values, in the limit n — > 1, of those of interest for us now. We have: 



d+(n) 



d 



and, defining 7y(n) ~ fij/{n — 1), 



4(n- 1) 



di, d_(n) 



16n/ 



w 48 w 



S(n)c 



27/V 

11 + h n s 

Po 

rif(n — 1) 



n f (n- 1) 



4n/(n-l) 
81 



(4.4a) 



Finally, 



r 12 = -l£n f T F C A , 
r 2 2 = ^fn f T F C A ~fn f T F C F . 



C[]\n) ~ C F [^ 2 - ^ - 4C(3)](n - 1), C&>(n) c fn/Tp 



C«(n) 



16T F (1) 127> 



3(n- 1) 



n- 1 



(4.4b) 



(4.4c) 



In the soft Pomeron hypothesis the behaviour of F$,g as x ^ is, as discussed, dominated by the singularities of the 
C^\n), 7^(71) as n — * 1, which in turn give those of the n{n). From (4.2 - 4) one easily finds, 



/ 



n{n,Q ) ~ 61a; 

n— >1 



-d+(n) 



1 + 



fe = 127> + 



2/3o 



n2 + r 2 2 + 



V n/(n- 1) 

A) 



aki 

n-l 

1 + - 



ak 

n- i/y 

3r i2 



(4.4d) 



, k\ = k — 



8n f 



Note that, because for n — » 1, d + (n) 3> d-(n), only the term in a s + (and not that in a s ) contributes. 
Next we evaluate 6 X in terms of Fj(x, Qq), assumed to behave as in (4.1) so that 



/•»»(«, Qo) - n c * 



n->l n — 1 



(4.5) 



This is accomplished using again (4.2) for Q = Q , profiting from the independence of b\ on Q 2 . The computation is 
straightforward and we find 



Air(n - 1) 



(4.6) 



7 Defined as in ref. 6, which fixes the arbitrariness in Fg- 
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Plugging the result into (4.4) we get 



1 + 
9 



3r i2 



8n/ 



H{n, Q 2 ) ^^coe Tda/i(n - 1) - Tdl 

\ 



n->l 

a s (Q 2 ) k _a s (Ql) 



47r(n - 1) 47r(n - 1) 



\n/(n-l) 



a s (Q 2 



<Xs{Ql) 



(4.7) 



'4n(n - 1) 47r(n - 1) J / 



r = log 

with do, di as in (4.4a). We then invert the Mellin transform. Generally, if 



and 



then 



where 8 



H(n) 



dx x n -"F{x), 



a d T/4(n-l) 



(n- iy 

x — >U 

1 >OC 

4 CT 

£ = logo;" , a = \v - §, / CT (r) = — (d T)~ 

7T2 



(4.8a) 

(4.8b) 
(4.9) 

(4.10) 



The proof is elementary and is obtained by substituting (4.9) in (4.8a) and integrating. Thus we get the final result, 

F s (x,Q 2 ) 



^r 1 



\64tt 2 J 



exp 



x 1 + 2 



d T 



ki 



a s (Q 2 ) 



47T 



k(l + 5) 



V doT^ - diT 

*s{Ql) " 



(4.11a), 



47T 



F G (x,Q 2 ) 



x a 



9c 



-exp 



\Mjt£ - diT 



d a r 



47T V ' 47T 



(4.11b) 



Here k, k\ are given in (4.4d) with the of (4.4b). We have added an arbitrary factor (1 + 5) for reasons that will 
be clear later. In the soft Pomeron model, one of course has 6 = 0. Numerically, 



h=k- 3r 12 /8n f ~ 42.19 
k ~ 22.19 (both for n f = 4). 



(4.11c) 



Eq. (4.11a) may be compared with the calculation of Ball and Forte. I 11 ! We agree in the LO term, and in the 
coefficient of a s (Q 2 )/47r in the NLO term, but disagree in the coefficient of the a s (QQ)/4:ir term. This is not of great 
moment 9 since the numerical difference is slight, 22.19 vs 16.19. Eq. (4.11b) is given here for the first time. 

NLO corrections are very large. Indeed, for fixed Q 2 , x — > 0, the NLO correction overwhelms the LO part. 
This, together with the problem posed by the BFKL-Florentine terms! 19 ' 20 ! 

i^ a ' ( " 2) ; uj = Constant 



8 From (4.8 to 10) it thus follows that powers of n — 1 correspond to powers of yj\ log:r|. Therefore, we cannot, unless 
a more definite assumption is made about the behaviour of the structure functions at Qq, give results more precise 
than terms of relative order 1/ yj log x \ . 

9 Nevertheless an independent calculation that resolved the discrepancy would be welcome. 
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will be discussed in Sect. 5. 



4.2. Theoretical calculations: longitudinal structure function. 



We define as before (Sect. 2.3) R' ~ F^/ Fs- Because, in the soft Pomeron dominated model, the contribution 
to Fl of Fs is subleading with respect to that of F G in the x — > limit, it follows that we may, to errors of relative 
size \/l/| logx|, neglect the contribution of Fs to Fl- For completeness, however, we will give the formula including 
this contribution of Fs- We then have an equation similar to (2.13), 



Jo 



F L (x,Q 2 )~F G (x,Q 2 ) dyC^(y,Q 2 )+F s (x,Q 2 ) dyC^(y,Q 2 ) 



/' 

Jo 



and the C L are as in (2.12b). 
We let 



I g (x)=T F n f [C F I G1 (0) + C A I G2 (0;x)] ~ -17.62 - 32n f^ CA log x ~\ 
I q (x) = C F {C A - 2C F )I S i(0) + C 2 F Is2(0) + C F T F n f I S3 (0) + C F T F n f I PS {Q;x), 



where 



dyc% )L (y) = n f T F [C F I G1 (0) + C A I G2 (0; x)}, 
and similar expressions for the Is, Ips- We then obtain the NLO expression for R' , 

R'(x,Q 2 ) ~ \2C F + l 



*0 



a s (Q 2 ) 
Air 



+48T F 



d T 



80 



T f Ca 



dor + * lg) 



a s (Q 2 )] \a s (Q 2 ) 



Att 



An 



(4-12) 



(4.13) 



(4.14) 



Because the NLO corrections are so large, we will use, instead of (4.14), a nonexpanded version for comparison 
with experiment. Removing also the NLO contribution of Fs the formula to be employed for numerical calculations is 
then, 



R'(x,Q 2 ) 



<5 2 >Q^ 



AST F \l-^ + 2C F 
dor 



1 



1 + 



{2k^/d^a s {Q 2 )/Air 
(2fci v ^7^) a s (Q 2 )/AiT 



(4.15) 



x\l + ll g (x) 



a s (Q 2 ) \ a s (Q 2 ) 



Att 



Att 



4.3. Comparison with experiment. 

For the soft Pomeron dominated model a very peculiar phenomenon occurs: the LO expressions produce fits 
better than the NLO ones. What is more, and unlike in the hard singularity case where we could blame the discrepancy 
on the large x points, here it is uniformly distributed. The strategy for comparison with experiment should be different 
now. First of all, we will not include a term like (1 — x) v connected with the saturation of the momentum sum rule 
since it is now very small, and would arrange nothing. Secondly, we give parameters for the LO fit for the restricted 
(x < 0.01) range, and we give results of the NLO calculation both for the restricted (x < 0.01) and full ranges. These 
we will discuss in greater detail. 

As stated we begin a LO calculation, fixing A(l loop Uf = 4) = 0.20 GeV, and taking for definiteness the 
HI data, plus the neutrino data for stability. We find the results of Table Vila. 
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Table Vila. LO calculation; x < 0.01, HI plus v data. 



(e 2 q )c (e 2 q )B NS Ql X>/d.o.i. 
0.12 0.44 0.46 GeV 2 



T 2 48.2 
68-3 



We consider next the NLO calculation which we split into two parts: restricted range, and full range. For the 
first we give the results of the calculation in Table VII&. Only the HI data are considered, for comparison with Table 
Vila. We do include neutrino data to force reasonable values for Bns- 



Table VII6. NLO calculation; x < 0.01, HI plus v data. 



A 



{e 2 )c 



(e 2 q )B NS 



Ql 



X 2 /d.o.f. 



0.160±^GeV 0.292±0.001 0.326 ± 0.020 0.86±^ 4 40 GeV 2 



74.3 

68-4 



As mentioned, the chi-squared has clearly deteriorated, 10 although the values of Qq, Bns are more realistic 
now. The value of A, also fitted, is reasonable. 



For the full range we give the results of the fits to both HI and Zeus data in Tables VIII. For the HI set there 
is no reasonable minimum for A; for Zeus the optimum is for A = 0.165 GeV. Thus we fix this value for both sets of 
data. 



Table Villa. NLO calculation; Zeus plus v data, x < 0.025, A = 0.165 GeV 



{e 2 q )c (e 2 q )B NS Ql x'/d.o.f. 



0.282 0.240 0.90 GeV 2 ^| 



Table VHIb. NLO calculation; HI plus v data, x < 0.032, A = 0.165 GeV 



(e 2 )c {e 2 q )B NS Ql x>/d.o.i. 



0.265 0.246 0.70 GeV 2 ^0Jl 



Our results are consistent among themselves, and with an existing NLO calculation, based on HI data^ 21 !; the 
comparison with the Zeus data for x < 0.025 is shown in Fig. 6. 



10 For Zeus data we would have had a much worse figure, x 2 /d.o.f.= ||^|, but a slightly better A = 0.20. The other 
parameters do not change substantially. 
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x = 0.00025 



10 12 14 16 rf- 



x = 0.0004 



10 15 20 25 n 2 



x = 0.00063 



10 18 26 34 



.r = 0.001 



10 24 38 52 



x= 0.0016 




10 34 58 82 



x = 0.0025 



a- 



10 50 90 130 n 2 



x = 0.004 



10 80 150 220 n 2 



x = 0.0063 




120 240 360 n 2 



x = om 




■jw-^H — 

4 





240 480 720 n 2 



x = 0.016 



330 660 990 



x = 0.025 



330 660 990 r,2 



Figure 6. Predictions for the structure function F% to NLO 
in the soft-Pomeron model, and Zeus data. 



The results may be confronted with the ones obtained if not including the NLO correction: we would have 
obtained, for the generally preferred value A = 0.23 GeV, 



X 2 /d.oi. = 1 ^ 2 ' 1 - (Zeus) 



X 2 /d-o.f. = 



120-3 
90.0 
110-3 



(HI). 



The situation is somewhat unpleasant. To make it worse, we mention that, if we delete the term in a s (Q§) in Eq. 
(4.11a), by simply putting S = — 1 there, the quality of the fit improves substantially: to a chi-squared/d.o.f. of 1 ^ Q 4 
for Zeus data, and for HI, with A = 0.23 GeV. 

It is difficult to draw a clear-cut conclusion from this. At any rate, in all cases the fits are comparable in quality 
to those obtained with the hard singularity hypothesis, and reasonably good; more discussion will be given in Sect. 5. 
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For the longitudinal function, the predictions 
and comparison with experiment are depicted in Fig. 
7; both the LO prediction based on the parameters of 
Table VII, and the NLO ones using the figures from 
Table Villa. Like in the hard singularity case, and for 
the same reason (large size of NLO corrections to Fs) 
there is a dramatic decrease between LO and NLO 
predictions, particularly for "large" values of x, and 
for very small ones. NLO results, depicted for various 
values of x in the figure are below the data. One 
cannot, nevertheless, consider the disagreement with 
experiment to be serious given the errors both of it 
and of the theory. Perhaps more serious is the problem 
that the NLO corrections also here overwhelm the LO 
piece for x — > 0. 

5. HARD PLUS SOFT SINGULARITIES. LARGE Q 2 . 
5.1. Discussion. 

It may appear strange that two mutually contradictory hypotheses, leading to so apparently different behaviours 
as the soft and hard Pomeron ones, produce both results in fair agreement with the data. The reason, however, is 
not difficult to find: both behaviours solve the QCD evolution equations, so the agreement of the calculations with 
experiment only depends essentially on the theoretical formulas fitting experiment at one value Q\ , say Q\ = 12 GeV 2 , 
and on the validity of QCD for the subsequent evolution for larger Q 2 at small x. However, neither the hard nor the 
soft Pomeron solution are fully satisfactory. The hard Pomeron expression fails to fit data with x > 0.01. The soft 
Pomeron does not produce a marvellous fit either, and in addition presents conceptual problems, that we now briefly 
discuss. 

First of all, we have the problem that, in the soft Pomeron case, the NLO overwhelms the LO term for small 
x (as <~ a s ^\ogx), so the soft Pomeron-inspired formulas must necessarily fail for the strict x — > limit. This is not 
the case for the hard singularity behaviour. Secondly, a power behaviour seems to be indicated for consistency with 
the 7*7 scattering case, where it has been shown to occur! 22 '. Finally, we have the following argument. By Reggeon 
calculus methods or summing ladder graphs, the authors of refs. 19, 20 have, at tree level, found the behaviour 

F 3 ~x-"-<»'\ U0 = 4C!a]o * 2 . 

This poses no threat to the hard singularity behaviour since this term is subdominant with respect to x~ x if, as seems 
natural, the argument of a s is proportional to Q 2 ; but it is incompatible with the soft Pomeron hypothesis because the 
new term dominates it. (However, there is a way out if the argument of a s was the hadronic energy, s ~ Q 2 /x for 
then, as x — > 0, the new term would become merely a constant). 

A possibility that allows us to keep the best of both worlds is that one has, at a low Qq, a behaviour sum of 
the hard and soft Pomeron ones, 

F s (Q 2 ,x)~c + bx- x . 

Although this implies that in the limit x — > the hard singularity will dominate, for finite x, if c <C b, both soft and 
hard singularities may contribute comparable amounts. In this context it may be remarked that a behaviour like the 
one above has been shown! 23, i] to describe very well photoproduction (Q 2 — 0) with constants precisely in the relation 
c<C i. If this persists up to Q 2 ~ a few GeV 2 , the mixed behaviour would be indicated. Needless to say, since 
both soft and hard-singularity dominated behaviours fit the data a mixed one will do so even better: for example, the 
deficiencies of the hard singularity picture at "large" x, and of the soft Pomeron one at all x, discussed in connection 



R(0,Q 2 ) 




80 GeV 2 160 GeV 2 240 GeV 2 



Figure 7. Predictions for R'(x,Q 2 ) to LO 
and NLO, soft Pomeron model. R'\ dotted, 
x — 0.05; dashed line: x = 10~ 3 ; continu- 
ous line: x = 10~ 4 . Hatched box: ref. 17 
(x ~ 10~ 3 ). Dots: data from ref. 18 at 
x ~ 0.05. 
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with the fit to the Zeus data, would likely be at least partially cured. 11 The conceptual difficulties of the soft Pomeron 
term alone also disappear as, in the strict limit x — > 0, it is dominated by the hard piece. 

There exists also a theoretical argument in favour of the hard plus soft Pomeron situation, and it comes from 
multi-Pomeron exchange theory. If, at a fixed Q 2 of the order of the GeV a single hard Pomeron gives 12 



F 1P (x,Q 2 ) ~ b lP {Q z )s x , 
s = Q 2 /x, 

then an n-Pomeron term will produce the behaviour 

F nP (x,Q 2 ) ~ b nP (Q 2 )s nX ; 

s — >oo 

the constants b nP should depend on the momentum at which they are calculated. 

In some approximations (e.g., of eikonal type! 24 !) one has, for Q 2 = —M 2 laA , i.e., for on-shell scattering of 

hadrons, 

b nP (-M 2 ad ) = (-ir +i -c, 

so we get for the sum 

oo 

F s (x, -M 2 ad ) = F nP (x, Q 2 = -M 2 ad ) 

n=l 

= C-C cxp[- KS A l ~ C. 

s^oo 

For values of Q 2 of the order of Qq ~ 2 — 4GcV 2 , we expect that the b nP (Q 2 ) will not change much from 
&np(-Mj 2 ad ), so if we write 

b nP (Ql) - bM-MLd) + A n, 

we will then get, 

F s (x, Ql) ~ F s (x, -M 2 ad ) + A n - C + A lX - x + A 2 x~ 2X + ..., (5.1) 

with the A small. This expression contains a hard plus a soft Pomeron of the type discussed above. It also contains a 
term A 2 x~ 2X , whose inclusion we will consider in Sect. 5.3. 

5.2. Hard plus soft singularities for F 2 . 

We next discuss the large Q 2 > 10 GeV 2 region, under various hypotheses for the small Q 2 region, which we 
then evolve with QCD. We will consider moderately large values of x, x < 0.032 because we will be interested not only 
on the leading behaviour as x — > 0, given almost certainly by a hard singularity, but also on the subleading corrections. 
So we assume that at a certain, fixed Qq ~ 1 GeV 2 , one has 

F s (x,Ql) ~ B s x- X +F corr .(x,Q 2 ). (5.2) 

For the correction term, -F C orr.(^, Qo), we consider the following possibilities: a soft Pomeron, 

^corr.( x 7 Qo) — constant; (5.3a) 

and a P' Regge pole, 13 

F^ IT {x,Ql) ~ constant x x 1 "^'^, a P -(0) - 0.5. (5.3b) 

This last possibility is considered because, as shown in ref. 10, any behaviour a > at Qq, produces at larger Q 2 
behaviours differing from the soft Pomeron one, (1.4a), only in the pre-factor, but with the same exponent. 

Once assumed the behaviours given in (5.3), and taking for simplicity that the gluon structure function behaves 
like the quark singlet one, we evolve with QCD for higher Q 2 . From the results of the previous Sects, we find, to NLO, 

F s (x, Q 2 ) ~ B S {l + CS(1 + J" )aS } c^ 1+x ^^[a s (Q 2 )]- d ^ +x ^x- x " + F COII .(x, Q 2 ), (5.4) 



11 The hard singularity picture overshoots the large x, Q 2 data (Fig. 3) while the soft Pomeron one undershoots the 

small x points, and undershoots large x ones: see Fig. 6. 
12 The following discussion is rather sketchy; details and references may be found in the review of ref. 24. 
13 For Regge pole theory cf. ref. 25. 
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and depending on the low Ql hypothesis we make we find the following correction terms: 



F^ r \x,Q 2 )^ 1 + 2 



9£t 



a s (Q 2 ) k a s (Q 2 y 



Air 



Air 



cxp [y/do^r- rfirj 



(5.5a) 



_4tt 2 (33- 2n/)_ 

for a soft Pomeron, Eq. (5.3a). For the P' Regge pole we find a very similar formula: 



^*corr.( a; > Q ) 



1 + 2 



h a s (Q 2 ) k a s {Ql) 



47T 



47T 



(5.5b) 



As we will see, none of the three possibilities gives a really good fit in the "large" 0.01 < x < 0.032 region; 
for the more precise Zeus data the xVd.o.f. is of 1.7. To remedy this we consider the possibility of softening the large 
x region by multiplying Fs{x,Q 2 ) by a factor (1 — x) v , as discussed in Sect. 3. Here, however, we take v constant 
because this already produces an excellent fit. 

The results are summarized in Table IX, where for definiteness we compare the fits obtained with (5.5a,b) with 
the fits found using only the soft Pomeron-dominated expression. We have not fitted A, which we have set equal to 
0.470, for reasons that will be apparent in next section; if we had fitted it, we would have obtained A = 0.43 and an 
improvement of only two units in the chi-squared 14 . We give the results for the Zeus data only; later we will present 
simultaneous fits to HI and Zeus data. 



Table IX.- n f = 4; Zeus data, plus neutrino data, x < 0.025, Q 2 > 12.5 GeV 2 



Soft Pomeron only 



Qo (e 2 )co 
0.165 GeV 0.90 GeV 2 0.282 



A 



(e 2 q )B NS xVd-o.f. 

0.240 T^7T~T 



* |A(f 



(fixed) 
47 



Hard + Soft Pomeron 

Hard + P'*, ( A (fixed) 
large x softened** | 0.47 

Hard + Soft*, ( A (fixed) 
large x softened 1 0.47 



Ql 

2.45 GeV 2 



(e 2 )c (e 2 q )B s 
0.252 4.42 x 10" 



(e 2 )B NS xVd-o.f. 
0.294 



Ql 
1.11 GeV 2 

Ql 

2.72 GeV 2 



{e 2 q )c P , 



{4)Bs 



0.616 4.25 x 10" 



(e 2 )B NS xVd-o.f. 
0.343 -r^-r 



(e 2 )c 
0.310 



(e 2 q )B s 
3.417 x 10" 



(e 2 q )B 



NS 

0.370 



X 2 /d.o.f. 

143 



* A fixed at 0.230 GeV. The optimum value would correspond to A ~ 0.45 GeV. 

** If we had not corrected for the large x values, i.e., we had not included the factor (1 

have obtained a xVd.o.f. of 270. 



x) v , we would 



In this table the expression "large x softened" means that we have multiplied the formulas for Fs by a factor 
(1 — x) u , v ~ 11, to correct the structure functions for (relatively) large values of x. For the hard singularity case, cf. 
ref. 5 and Sect. 3.4 here; for the Hard + Soft singularities case, we have taken v = 10 ~ 11. (We will discuss further 
the "large" x region in Sect. 5.3). 

5.3. Hard plus soft singularities, plus triple Pomeron term for F 2 . Best (global) fits. 



To be precise, if we had fitted A to e.g., the Zeus data using a hard plus soft term, we would have obtained (not 
correcting for large x, and fixing A = 0.23 GeV), 

A = 0.429, Qo = 2.40 GeV 2 , (e 2 )c = 0.244, {e 2 )B s = 3.83 x 10~ 4 , {e 2 )B NS = 0.30, 

for a X 2 /d.o.f. = T ii^. 
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Clearly, the best fit is obtained with the hard plus soft Pomerons. Not only the x 2 /d.o.f. is quite good, but 
the values of the parameters are very reasonable. In fact, more evidence in favour of the "hard plus soft" scenario will 
be given in next section; for now we will consider that the fits are so good, that it makes sense to see if one can find 
evidence for a "triple Pomeron" term. That is, we consider that [cf. Eq. (5.1)] 

F s = F Soft + F H ard + F TP , 



F TP (x,Ql) ~ (Const.) x 
so that, when evolved with QCD to large Q 2 , 

c s (l + 2\)a 



-2\ 



F^{x,Q 2 )^B TP ll + 



47T 



-2A 



(5.6) 



Note that this is 0(af +< - 1+A ' ) d +( 1+2A )) ~ 0{a], A ), i.e., subleading in powers of a s , with respect to Fn a rd- 

We present in Table X the parameters of the fits to, simultaneously, HI and Zeus data on ep, plus neutrino 
data. This gives our best set of formulas, providing an excellent fit to experiment in a very wide range of Q 2 , x. In 
the second case (Table X6) we do not give the fit including a triple Pomeron term as the x 2 /d.o.f. does not vary 
appreciably if including it provided (c 2 )\Btp\ ^2 x 1CP 4 . We consider the parameters given in Table Xa and Table 
Xc (see below) to be the more reliable ones for describing small x structure functions. If we had fitted also A with 
the whole set of data we would have obtained minima for values comprised between 0.42 and 0.49, with a variation of 
the chi-squared of less than two units with respect to the one obtained fixing A = 0.470. Finally, if we fit the QCD 
parameter A, the values which provide minima vary between 0.555 GcV and 0.310 GcV, and the chi-squared improves 
by less than five units. Because of this we consider, as stated, that it is justified to favour the fits obtained with fixed 
A = 0.470, A = 0.23 GeV. 



Table Xa. n } = 4; Zeus plus HI data; Q 2 > 10 GeV 2 , x < 0.01. 



Hard + P 



{ A 2 



(fixed) 
,47 



Ql 
2.95 GeV 2 



(e 2 )c P 



{4) B s 



Hard + P, 
+TP term 



A (fixed) 
0.47 



Ql 



4.45 GeV 5 



<e 2 )c P 
0.258 



0.296 4.28 x 10 
8.33 x W~ 4 



{e 2 q )B NS 
0.349 



{e 2 q )B T p 
1.67 x 10" 



xVd-o.f. 

138.3 
144-4 

{e 2 q )B NS ; 
0.359 



K 2 /d.o.f. 

129.3 
144-5 



Table Xo. n f 



4; Zeus plus HI data; Q 2 > 10 GeV 2 , x < 0.032. 



Hard + P 
x "softened" with ( 



fA(f 
1 -x) v \ 



(fixed) 
47 



Ql 

2.28 GeV 2 



(e 2 )cp 
0.311 



(4)Bs 
2.72 x 10" 



(e 2 q )B NS xVd-o.f. 
0.315 



Table Xc. n f = 4; Zeus plus HI data; Q 2 > 10 GeV 2 , x < 0.032.. 



Hard + P 
+P' 



A (fixed) 
0.47 



Ql 
5.00 GeV 2 



0.588 



(e 2 q )C P , 
-0.271 



{e 2 )B s 
4.66 x 10" 



(e 2 )B NS xVd-o.f. 
0.262 ^4 



A fixed at 0.230 GeV. NLO corrections included 



We discuss now in some detail the larger x region. If taken by themselves, both soft and hard Pomeron 
expressions (and a fortiori a sum of the two) must, as discussed in ref. 1 and Sect. 3.1 here, run in contradiction with 
the momentum sum rule if mantained for fixed x and Q 2 — > oo; and this contradiction starts becoming noticeable at 
the higher Q 2 ~ 1 000 GeV 2 for x > 0.02: so a modification of our formulas for finite x is necessary. In the present 
paper we have, until now, introduced it phenomenologically by multiplying the low x expressions by a factor (1 — x) v 
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( "softening"), with v constant or depending on Q 2 . A more rigorous procedure would be to assume, at a fixed Qq, a 



behaviour like 



F s {x,Ql) = (B s x- x + C){l-xY\ 
and then evolve with QCD. This is best done by expanding first (5.7) in x, 

F s (x, Ql) ~ B s x~ x +C- v Q B s x 1 - x + .... 



(5.7) 



(5.8) 



The dots correspond to terms behaving, as x — > 0, as higher powers of x, which need not be considered; see below. 
(5.8) has exactly the form of a hard Pomeron, plus a soft Pomeron, plus a P' Regge pole, plus higher powers of x. 
From the results of the previous sections we know that all terms vanishing for small x yield the same expression, up to 
a constant, when evolved to large Q 2 , as the P' piece. So we may lump the piece —vnBsx 1 ~ x + . . . into a single 
term like that of Eq. (5.5b), to be added to a soft and a hard term. To LO thus, 



F s (x, Q 2 ) ^ B s a; d + + F c p OII ix, Q 2 ) + F%„Xx, Q 2 ) 



(5.9) 



F^ TT (x,Q 2 ), F^ orr (x,Q 2 ) given in Eqs. (5.5). The resulting fit (including NLO corrections) is described in Table Xc 
above. This certainly improves the fit at large x with respect to the unsoftened situation, but not much; indeed, less 
than the simple "softening" used before in the text 15 , and it certainly does not solve the momentum sum rule problem, 
either. This should not be too surprising: by its very nature, a calculation with leading terms only for x must fail for 
larger values of this variable. 

What we mean by this exactly is the following. Consider that at a fixed Qq one had exactly a hard Pomeron, 
Fs(x,Qq) = Bx~ x . Then, at any larger Q 2 , we have the moments [in this simplified discussion we neglect the matrix 
character of the evolution equations] 



<*.(Q 2 ) 



a s (Q 2 ) 



d + (n) 



n(n,Q 2 ) = B 

Writing identically c Td +^ = C r<i +( 1 + A )[l + 5{n, A)] we find 

F s (x,Q 2 ) = F^ -rr s , r s 
and the subleading piece is such that it has moments 

V SL (n,Q 2 ) = B s {a s (Q 2 T d +^ 



(1 + A)' 



2 , _ pLcad. + F SL. F Lcad. (a , ; g 2) = Bg ^ (Q 2 -d+ (1 + X) x ~ A ^ 



6(n, A) 
n- (1 + A)' 



S(n, A) vanishes for n = 1 + A; hence the first singularity of the ^ SL (n, Q 2 ) occurs for n = 1 
and there, 



V SL (n,Q 2 ) ~, 



-B S 



, K d + (i+\) expr 



♦i 1 + A 

This is of the soft-Pomeron type apart from the factor a s 

\x,Q 2 



a -d + (l + X)+d le 




If we continie to subtract the terns singular at n = 1, 2, . . ., we would get an asymptotic series for Fg- For x — > this 
is dominated by the term F^ cad -- but for fixed x, the remaining terms in the series end up by overwhelming i^ead. as 
Q 2 becomes very large. 

A precise evaluation is not difficult; it would also involve the gluon component. We shall not present the 
corresponding fits here; to do so one would have to include also subdominant corrections to the soft Pomeron piece. It 
is unclear that the effort would be worth the results, given the good quality of the fits with dominant terms only. 



'Nevertheless, from the point of view of rigorous QCD the softening given in Eq. (5.9) is to be preferred to mere 
multiplication by (1 — x) v for all Q 2 ; e.g., for extrapolations to higher Q 2 , since (5.9) is compatible with QCD 
evolution for x — » 0. What is more, the parameters are fairly stable from Table Xa ot Table Xc, advantages that 
in our opinion offset a small increase in xVd.o.f. . wit respect to Table X6. 
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Figure 8a. Comparison of predictions from Table X with HI ep data for F2. 
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Figure 8b. Comparison of predictions from Table X with Zeus ep data for F2. 



- 23 - 



-Small x deep inelastic scattering- 



In Figs. 8 we show the comparison of our fits, with the parameters in Table X, with data. Note that the same 
values of the parameters are used in Fig. 8a and Fig. 8b. For both Figs. 8 we give the fits with the "softened" and 
straight formulas: the continuous lines indicate large- x softening, and the dotted lines no softening. The improvement 
in the quality of the fits when compared with the "soft-Pomeron only" or "hard-Pomeron only" is obvious to the naked 
eye here. 

5.3. Gluon and longitudinal structure functions. 

Detailed predictions for the gluon and longitudinal structure functions are obtained trivially by adding the soft 
and hard Pomeron expressions given in the previous sections We will leave the details of this to the reader; likewise, we 
do not draw the figures for the Fq, Fl as they would not differ much from the ones drawn for, say, the hard Pomeron- 
only hypothesis. However we would like to comment here on a particularly interesting prediction of our analysis for the 
growth of the cross-section <t 7P ^j/^ p (VK) as a function of the cm. energy, W. In fact, this cross section may be 
expressed as a function of the gluon structure function F G , 

v^j/toiW) = AFq(x = a-0f, Q 2 = M 2 ^) ~ 9GeV 2 , 

and A, a are constants approximately known. For Fq we have, 

F G (x,Q 2 ) ~ B G [a s ]- d +x- x 

and Bq may be calculated in terms of A, B$- So, using our formulas we have, for the logarithmic slope of the cross 
section, 

= log o-^j/^p (HQ ^ 2A 

logVF s^oo 

The figure reported in a fit! 26 ! including recent HERA data^ 27 ! gives 5 = 2A = 0.9. This in very good agreement with 
the optimum values of A obtained with the hard plus soft fits, 2A = 0.83 to 1.0, and is clearly superior to the results 
following from the hard Pomeron only hypothesis, 2A = 0.64 to 0.76. 

6. HARD PLUS SOFT SINGULARITIES: SMALL Q 2 . 



The quality of the results obtained by assuming that at values of Qq ~ 3 GcV 2 one has a hard singularity, 
x~ x , plus a soft (constant) Pomeron term, evolved with QCD to large values Q 2 > 10 GcV 2 , leads us naturally to 
the question wether it is possible to extend the analysis to the low Q 2 region as well, thus enabling us to address the 
important issue of the connection between the perturbative regime (Q 2 > 10 GeV 2 , say) and the region Q 2 < 10GeV 2 
where nonperturbative effects are determinant. 

It should be obvious that, unless one were able to perform a full, nonperturbative calculation, we must content 
ourselves with phenomenological evaluations. Here we use approximate, QCD-inspired formulas and assumptions and 
enquire wether we can still fit the data. We will find that this is indeed the case; in particular, we will see that the 
extension of the fit of the data to Q 2 — > implies self-consistency conditions both for the singlet and the nonsinglet 
which will allow us to calculate the constants A, p, getting values in impressive agreement with other (in particular, 
high Q 2 ) determinations. 

The expression for the virtual photon scattering cross section in terms of the structure function F 2 is 

er 7 (Q2 =0)p (s) = ^-F 2 (x, Q 2 ), with s = Q 2 /x. (6.1) 

We would like to describe this down to Q 2 — ► 0. In the low energy region we should, as discussed, take the soft-Pomeron 
dominated expression to be given by an ordinary Pomeron, i.e. , behaving as a constant for x — > (or equivalently, 
s — > oo): the expression for F 2 that will, when evolved to large Q 2 yield (5.4), (5.5a) is 

F 2 = (c 2 q )\B s [a s (Q 2 )}- d ^ 1+x ^x-^ 

, (6-2) 
+C + B NS [a s (Q 2 )]- dNsil - pa) x pa j. 
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Because NLO corrections are large for Q 2 < 10 GcV 2 and we are interested in a semi-phenomenological description, 
only LO formulas will be used. Note also that the C in Eq. (6.2) is different from the Co in, say, (1.4a) as the gluon 
component also intervenes in the evolution. 

On comparing (6.1) and (6.2) we see that, as noted in ref. 1, we have problems if we want to extend (6.2) to 
very small Q 2 . First of all, 

- w> > " (6 - 3) 

diverges when Q 2 ~ A 2 . Secondly, Eq. (6.1) contains the factor Q 2 in the denominator so the cross section blows up 
as Q 2 — > unless F 2 were to develop a zero there. 

It turns out that there is a simple way to solve both difficulties at the same time. It has been conjectured I 28 l 
that the expression (6.3) for a s should be modified for values of Q 2 near A 2 in such a way that it saturates, producing 
in particular a finite value for Q 2 ~ A 2 . To be precise, one alters (6.3) according to 

as[Q ] l3 log(Q 2 + M 2 )/A 2 ' 

where M is a typical hadronic mass, M ~ m p ~ A(rif = 2) . ..; the value M — 0.96 GeV has been suggested on 
the basis of lattice calculations. It has been argued that saturation incorporates important nonperturbative effects. In 
the present paper we will simply set M = A = A cS , to avoid a proliferation of parameters. For the Pomeron term 
[the constant in Eq. (6.2)] we merely replace C — ► Q 2 /(Q 2 + A 2 S ), using a procedure similar to that of ref. 29. The 
expression we will use for low Q 2 is thus, 

F 2 = (e 2 ){B s [a s (Q 2 )]- d ^ 1+x '>Q- 



2A S A 



+C Q2 + A* + BNs[a s {Q 2 )]- d » s{1 - p) Q 2p s-e), 



(6.4a) 



where 



47T 

) = Tl ira i A2 ma2 ( 6 - 4b ) 
l3 Q log{Q 2 + A 2 cS )/A 2 cS 

and we have changed variables, (Q 2 , x) — > (Q 2 , s = Q 2 /x). 

We have still not solved our problems: given Eq. (6.1) it is clear that a finite cross section for Q 2 — > will 
only be obtained if the powers of Q 2 match exactly. This is automatic by construction for the Pomeron term, but 
for the hard singlet and the nonsinglet piece it will only occur if we have consistency conditions satisfied. With the 
expression given in (6.4b) for a s it diverges as Const. /Q 2 when Q 2 — > 0: so we only get a matching of zeros and 
divergences for ct 7 (q2 =0 ) p (s) if A = An, P = po sucn tnat 

d+(l + A ) = 1 + An, d NS (l -p ) = l- Pa . (6.5) 

The solution to these expressions depends very little on the number of flavours; for nf = 2, probably the best choice at 
the values of Q 2 we will be working with, one finds Ao = 0.470, po = 0.522. The second is in uncanny agreement with 
the value obtained with either a Regge analysis in hadron scattering processes, or by fitting structure functions in DIS. 
The first is larger than the value obtained in the fits to DIS with only a hard Pomeron, which gave A = 0.32 to 0.38; 
but falls in the right ballpark of values obtained in the previous section with hard plus soft Pomeron, A = 0.43 toO. 5. 
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Figure 9. Comparison of predictions with data, Zeus plus HI. The neutrino data and 
prediction are not shown. We plot F2 vs. x 1 ^ 2 . 



We are perfectly aware that, by using Eqs. (6.4) down to Q 2 = 0.32 GeV 2 we are pushing perturbative QCD 
well below its region of applicability, and that the condition of matching at Q 2 — > is at best only of phenomenological 
value. Nevertheless, the fact that we get such reasonable predictions for Ao, po probably indicates that our procedure 
represents, grosso modo, the actual situation, which is also justified by the quality of the fit Eq. (6.4) provides. If we 
take all HI and Zeus data for 0.31 GeV 2 < Q 2 < 8.5 GeV 2 , and we include also 10 neutrino xF 3 data we find 

A cB = 0.87 GeV, (e 2 q )B s = 5.28 x 10~ 3 , (e 2 q )B NS = 0.498, (e 2 q )C = 0.486, 

for a x 2 /d.o.f. = j§§^- The value of A c g we have obtained lies somewhere inside the expected bracket, A(nj = 2) ~ 
0.35 GeV and the value found in the quoted lattice calculation for M, 0.96 GeV. Clearly, the fit gives a compromise, 
phenomenological quantity. 

The agreement between phenomenology and experiment, shown graphically in Fig. 9, is unlikely to be trivial; 
x varies between 6 x 10~ 6 and 4 x 10~ 2 , and F2 changes by almost one order of magnitude. To see more clearly this 
nontriviality, we replace the hard singularity by an evolved soft Pomeron, with a saturated a a . That is, we now fit with 
the expression 
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F* = {%)> 



Co 



^log[ as (Q 2 Q )/a s (Q 2 )] 



2\n 4 



4tt 2 (33- 2n f ) 



+C 



cxp ^ 

Q 2 



'dot 



S 5 S (Q 2 ) 



Q 2 + K s 



+ B NS [& S (Q 2 



iNs(l-p)Q2p g -p I 



a, as before. Then we find 



A cS = 0.41 GeV, (e 2 q )c Q = 0.094, (e 2 q )C = 0.253, (e 2 q )B NS = 0.41 

and a much deteriorated x 2 /d.o.f. = 250/(104— 4). We consider these results as convincing proof of the necessity of 
a hard component also at low Q 2 . 



7. DISCUSSION 



The main outcome of our analysis in the present paper is that we are able to give a unified, consistent 
description of small x DIS data, both for large and small values of Q 2 , by assuming, at low momenta, the presence 
of a hard plus a soft Pomeron, a procedure which improves substantially the quality of the descriptions with only one 
of these. There also appears some evidence for a triple Pomeron contribution; evidence which is, however, somewhat 
marginal. Besides this, there are a number of specific points to which we would like to draw also attention. 

First of all there is the matter of the dependence of our low Q 2 results on the saturation hypothesis for a s . It 
is clear that the good quality of the fits indicates that, with suitable modifications, perturbative QCD supplemented by 
saturation, may give a phenomenological description of the data down to very low momenta; but of course this should 
not be construed as a proof of saturation, in particular of the very specific form considered here. One may interpret 
our results, however, as showing that the saturation expression is particularly adapted to represent, in DIS, a variety of 
effects: higher twists, renormalons, and likely also genuine saturation. 

A second question is the connection between low and high Q 2 . For the hard piece there is no problem, as 
both expressions are identical up to NLO corrections. For the soft piece, if we start with a constant behaviour for 
Q 2 ~ 2 — 5 GeV 2 , then as Q 2 grows an expression like (2) will start to develop, the details of this will depend on 
what one assumes for the gluon structure function. Because the variation both with Q 2 and x of the soft piece is slower 
than that of the hard part, we think the best procedure is to assume constancy of the soft piece up to Q 2 = 8.5 GeV 2 , 
and the evolved form from there on; since a very good fit is obtained at the low momentum region already with the 
constant behaviour there is little point in adding frills, and a new constant (the soft component of the gluon structure 
function). 

Next, we say a few words on the parameters, starting with the QCD parameter, A. It is impossible to give a 
reliable determination of the value of this parameter from low x data alone; if fitting it, the central values vary between 
0.08 GeV and 0.55 GeV. If we take what we consider the more reliable fits, those in Tables Xa, c, and allow A to 
vary, we find an optimum value of 0.31 GeV; in particular, from the fit with hard plus soft Pomerons, plus P' , the 
optimum is A = 0.32 ± 0.05 McV for a x 2 /d.o.f. of 2 2 ^ 6 , hardly improving the result reported in Table Xc. This 
tells us little more than rough compatibility between the low x and other determinations of A. 

With respect to other parameters, we can say that the fits give precise determinations of the nonsinglet 
parameters, p, -B/vsl but the singlet parameters are much less precisely determined. For example, A varies from 0.38 
(LO, hard singularity only) to 0.32 (NLO, hard singularity only) to 0.47 ± 0.04 (hard plus soft, the best value in our 
opinion). Likewise, Bs varies by almost one order of magnitude. The reason for this may be traced to the dependence 
of the parameters on the theoretical formulas used to fit the data, in particular when going from LO to NLO because of 
the large size of the NLO corrections to Fs- We consider the parameters given in Tables Xa, c to be the more reliable 
ones in particular for extrapolations to larger Q 2 . 

APPENDIX 



In this Appendix we present the full collection of formulas necessary to evaluate electroproduction to NLO. 

Leading order quantities. 



-Small x deep inelastic scattering- 

/3 = ll-|n / , /3^102-fn/. 

d NS (n) 



^(«). D(n)= ^ (0) W 



2/?o 



Here, 7ws( n ) = 7nV) and 



(0) _ _32 
' — 3 



n 2 +n + 2 33 - 2n/ 9 
V 2n(n 2 - 1) 16 +I 



2/?o 



n 2 + n + 2 
n(n + l)(n + 2) 
1 1 



(A.l) 



We define, generally, the functions 



fe=i 



1 



k l (k + n) 1 



n(n-l) (n + l)(n + 2) 
S+(in) = 5 ; (n/2); S r (In) = S, 



5i(n) 



(A.2) 



n — 1 



and 



^(n) = -fC(3) T 71-^ Si{k + n). 



k=l 



The matrix that diagonalizes D(n) is S(n), 

S(n)-D(„)S(„) = (<♦<»> ,_»„,). 

with the eigenvalues ordered so that rf + > gL. S may be written as 

/ x gi2(n) 

ar„\- d-(n)-d+(n) 
W ~ rf + (n)-Dn(n ) d_(n)-£>n(n) 

V D 12 (n) 



(A.3) 



d-(n) — d+(n) 
Nonsinglet NLO quantities. 

We will not give explicit formulas for the Cj^g(n), J^ s {n), which may be found, misprint free, in refs. 6, 7. 
We will only present a few figures for relevant values of n from which good interpolation formulas may be written. One 
has 

C$1 (n = 1) = 0; 7 «(1) = f [13 + 8C(3) - 2n% 

the last formula for crossing-even functions (like Fjvs in electroproduction). For crossing odd functions, like xF 3 in 
neutrino scattering, the corresponding NLO anomalous dimension coefficient verifies 

7^,odd( 1 )=0- 

For n near 0.5, 







C^ s (n) = 


n = 0.4 
17.1; 


n = 0.5 
9.6; 


n = 0.6 
5.5 










rif = 3 




n f =A 




nf = 5 




n = 


0.4 


0.5 0.6; 


0.4 


0.5 


0.6; 


0.4 0.5 


0.6 


7& = 


-273.7 


-159.3 -98.2; 


-271.7 


-155.9 


-95.0; 


-269.6 -152.5 


-91.7 



Note that cff s (n) is independent of n/. 
Singlet NLO quantities. 

The four quantities Cij(n) may be found in ref. 6. Here we give only the two that enter the calculation for 
ep scattering. With C F = 4/3, C A — 3, Tp = 1/2, 



C[\\n) = C F ^2S 2 (n) + 3Si(n) - 2S* 2 



(n) 



2Si(n) 
rc(n + 1) 



« 3 4 
9 + - + 



n n + 1 n 2 J ' 



(A.4a) 
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,2 



dV(n) = AntTp [-- + !- + -£ - SUn) -^—^r^—\ ■ (A.4b) 

12 w ' \ n n 2 (n + l)(ra + 2) v ' n(n + l)(n + 2) J v 7 

Finally, for the NLO anomalous dimension matrix 7^ 1 ^(n) we have the following expressions: 

7^V) = C% |l65 l( n) n2 2 ( "+ 1 1)2 + 16 ( 2 5 l( n) - [5 a(n) - 

+24S 2 (n) + 645(n) - 8S 3 +(±n) - 3 - 8 3 "* + ~ 1 - 16 2n \. - 

3 2 n 3 (n + l) 3 n 3 (n+l) 3 

+CfCa {f^(n) - 8 (2Sx(n) - ^yy) [25 2 (n) - S+(|n)] 

-f 5 2 (n) - 32S(n) + 4S 3 +(±n) - f (A ' 5a) 



4 151n 4 + 236n 3 + 88n 2 + 3n + 18 2n 2 + 2n + 1 

+ 8 
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+n f T F C F {-if + fS 2 (n) + § 



9 

16 lln 7 + 49n 6 + 5n 5 - 329n 4 - 514n 3 - 350n 2 - 240n - 72 
+ ~ (n- l)n 3 (n+l) 3 (n + 2) 2 



-7#(n) = 8nfT F CA { [ - 2S 2 (n) + 2S 2 (n) - 2S+(±n)] 



n 2 + n + 2 
n(n + l)(n + 2) 
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^ L J n(n + l)(n + 2) n 2 

lln 4 + 26n 3 + 15n 2 + 8n + 4 1 
+ n 3 (n + l) 3 (n + 2) J ' 

This corrects a misprint in ref. 8 ( a figure 262n 3 instead of 26n 3 in the third line). 

■n + 2 
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-Small x deep inelastic scattering- 
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n(n-l) (n+l)(n + 2) 
148 32 32 



+ 



64 _ 

3 n-1 9(i(n+l) T n + 2 (n - l) 2 
256 32 32 64 1 



3n 2 3(n+l) 2 3(n + 2) 2 n 3 (n + 1) 3 (n + 2) 3 J 

This last corrects two misprints of ref. 8, a factor UfTp = Tr instead of Ta in the first line, and a sign, 
+l/(n + l) 2 instead of -l/(n + l) 2 in the fourth line. 

We finish by giving two values of the \ 7^ and tables with a few listings, sufficient for the calculations 
we are interested in. As for the first, we have 

C${2) = i C ( ${2) = -\n f , 



7 (1) (2) =2 1 



F s {x,Q 2 ) ~ B S {1 



1833 „ 



367 - 39n/ 
-(367- 39^) i|fn/ 
We give the tables for the quantities relevant for the exponential expression, qs and eg, 

c s (l + \)a s (Q 2 )- 

47T 
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0o 2/3 2/3 (d- - d+ + 1) d_ - d+ 
for the expanded expression just note that ws = qs + C5. Then, 
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For the gluon structure function, 



F G (x,Q 2 ) ~ B s 
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Coefficients and integrals for the longitudinal structure function. 

We give the coefficient functions for ep scattering, with unified notation. With the definitions of Eqs. (2.12), 
and with the polylogarithm functions 16 



we have, 
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Li,2 denoted £1,2 in refs. 12, 13. 
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-Small x deep inelastic scattering- 



<4 1)L (x) = 2n f T F C F | (8 + 24x - 32x 2 ) log(l - x) 
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Integrals: for n/ 


= 4, 
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(the values of the integrals not given when x-dependent). 
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